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Abstract 

In this paper, we study the polynomial stability of analytical solution and convergence of the 
semi-implicit Euler method for non-linear stochastic pantograph differential equations. Firstly, 
the sufficient conditions for solutions to grow at a polynomial rate in the sense of mean-square 
and almost surely are obtained. Secondly, the consistence and convergence of this method are 
proved. Furthermore, the orders of consistence (in the sense of average and mean-square) and 
convergence are given, respectively. 
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1 Introduction 

Stochastic pantograph differential equations(SPDEs) arise widely in control, biology, neu¬ 
ral network, and finance, etc. Asymptotic stability of analytical solution has received consider¬ 
able attention in literature for both deterministic and stochastic functional differential equations. 
Especially, plenty of literature on stability exist with non-exponential rates decay of the solu¬ 
tions (see [6] [9] [To] [Ej). One important non-exponential rates of decay is polynomially asymptotic 
stability, which means that the rate of decay is controlled by a polynomial function in mean-square 
or almost surely sense. This type of stability has been studied in m, [3i and |1|. Buckwar and Ap¬ 
pleby consider the polynomial stability of one dimensional linear stochastic pantograph differential 
equation in |2|, where the sufficient conditions of polynomially asymptotic property are given. 

The convergence of numerical method is another crucial property of stochastic differential 
equations. Recently, many researchers devoted to the stochastic delay differential equations. Mao 
Wei m gave the sufficient conditions of convergence with semi-implicit Euler method for variable 
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delay differential equations driven by Poison random jump measure. Fan studied the approximate 


solution of linear stochastic pantograph differential equations with Razumikhin technique in [8]. 
Baker and Buckwar, in [5], also investigated the linear stochastic pantograph differential equations, 


and the sufficient conditions of convergence of semi-implicit Euler method are obtained. 

This article investigates non-linear stochastic pantograph differential equations. It is organised 
as follows. In section 2, we introduce necessary notations and results on pantograph differential 
equations. The sufficient conditions of polynomial stability of analytical solutions in the sense of 
mean-square and almost surely are given in section 3. Finally, both the consistence and convergence 
of the semi-implicit Euler method are proved, and the orders are achieved. 

2 Preliminary notations and properties of pantograph differential 
equations 

Throughout this paper, unless otherwise specihed, the following notations are used. If x, y are 
vectors, the inner product of x,y is denoted by {x,y)=x'^y, and |x| denotes the Euclidean norm of 
X G We use to denote the transpose of A, if A is a vector or matrix. And if A is a matrix, 
its trace norm is denoted by |A| = y/traceA^A. 

Let (n, P, F) be a complete probability space with a filtration {Ft}t>Q satisfying the usual con¬ 


ditions. Let 0 < q < 1, x{t) = (xi(t), X 2 (t), • • • , X(i{t)Y' S and B{t) = {Bi{t), B 2 {t), • • • , Bm{t)Y' 
be a m-dimensional Brownian motion. Suppose h{x) be a function, the derivative of h{x) is dehned 
as D~^h{x), i.e 


h{x -|- 5) — h{x) 


D^h{x) = limsup 
5^-0 


( 2 . 1 ) 


Before studying the stochastic pantograph differential equations, the properties of deterministic 
pantograph differential equations are hrstly introduced. The equation has the following form: 


x'{t) = ax{t) + bx{qt) 
x{0) = xo- 


( 2 . 2 ) 


Lemma 2.1. Assume x{t) he the solution of (2.2), xq > 0, if a < 0, then there exists a constant 
number C such that 



(2.3) 


where a G R and a + bq°‘ = 0. 


Remark 2.1. If there exists a constant number C > 0, such that 


x{t)\ < C'|x(0)|t“, t >0, 


(2.4) 


then (2.3) holds. 
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Lemma 2.2. Assume x{t) be the solution of equation (2.2), xq > 0. If b > 0, p{t) : R~^ —)■ R~^ 
is a non-negative continuous function satisfying 

D^p{t) < ap{t) + bp{qt), t >0, (2.5) 

and 0 < p{0) < x(0), then p{t) < x{t) for all t > 0. Furthermore, if d < 0, and p{t) satisfys (2.5), 
then there exists a constant C > 0 such that 


where a + = 0. 


p{t) < Cp{0)t°‘,t > 0, 


( 2 . 6 ) 


3 Polynomial stability of analytical solutions for non-linear stochas¬ 
tic pantograph differential equations 


In this section, we consider the following equation 

f dx{t) = f{t,x{t),x{qt))dt-I g{t,x{t),x{qt))dB{t), 0 <t <T 
\ x(0) = XQ. 

where / : [0,T] x R^ x R^ ^ R^, g : [0,T] x R^ x R‘^ ^ /(t, 0,0) = 0, g(t,0,0) = 0, and 

E\xq\^ < oo. It is easy to see there exists zero solution for (13.71) . 


Definition 3.1. The zero solution for {3. ?[ ) is said to be mean-square polynomial stable, if there 
exists a constant number a < 0, such that 


lim sup 

£—)-oo 


logE\x{t)\‘^ 

logt 


< a, 


(3.8) 


where x{t) is the solution for |g. 7]) with any initial value x(0) = xq. 


Definition 3.2. The zero solution of (3.1) is said to he almost surely polynomial stable, if there 
exists a constant number a < 0, such that 

r log\x{t)\ 

iimsup—;- <a a.s., (3.9) 

logt - ’ ^ ^ 

where x{t) is the solution of {3. ?D with any initial value x(0) = xq. 


Assumption 3.1. Assume there exist real numbers a, b > 0,c > 0,d > 0, such that the coefficients 
f and g satisfy 


(xi - X 2 , fit, xi,y) - fit, X 2 ,y)) < a|xi - X 2 l^, Vxi, X 2 , 2 / £ R!^,t> 0; 


\fi't,x,yi) - /(t,x, 2 / 2)1 < b\yi - 2 / 2 I, 
\git,xi,y) - git,X 2 ,y)\ < c\xi -X 2 I, 
\git,x,yi) - git,x,y 2 )\ < d\yi - 2 / 2 !, 


Vx, 2/1,2/2 e R‘^,t> 0 ; 
Mxi,X2,y G R‘^,t > 0; 
Vx, 2/1,2/2 € R‘^,t> 0 . 
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According to /(t,0,0) = 0, 0,0) = 0 and assumption 13.91 we can estimate 

{x,f{t,x,y)) < {a + ^)\x\‘^ + ^\y\'^, \g{t,x,y)\‘^ < 2c‘^\x\'^ + 2(f\y\‘^. 

Theorem 3.1. Suppose that x{t) is the solution of ^3. ?[ ), if the coefficients f,g satisfy assumption 
roi and o + 6 + < 0, then the zero solution of is mean-square polynomial stable. 

Proof. According to definition 3.1, we just need to prove that there exist constant number 
Cl and a < 0, such that E\x{t)\’^ < CiC. 
ltd formula shows that 

= E{\xq\^) + E{ f [2{x{s), f{s,x{s),x{qs))) + | 5 (s, x(s), x(gs))p]ds). (3.10) 

Jo 

Let Y{t) = E{\x{t)\‘^), then for any t > 0, t + h > 0, we have 




1, 


1. 


rt-\-h 


Y{t + h) -Y{t) < 2E{ / [(a + -b)\x{s)f + -b\x{qs)f]ds) + 2E{ / [c \x{s)\ + d \x{qs)\ ]ds) 

Jt ^ ^ Jt 

rt+h rt+h 

<(2a + 6 + 2c^) / y(s)ds + (6 + 2d^) / y(gs)ds. 

Jt Jt 


Due to 


/o+'*y(s)ds- /Qy(s)ds 


Y (t) = limsup 
h^O 


Y (qt) = limsup 

h^O 


lY ^ ~ Jo ^ 


thus 


D+Y(t) <(2a + b + 2c^)Y(t) + (b + 2d^)Y(qt). 


(3.11) 


Note that 2a + 6 + 2c^ <0,6 + 2d^ > 0, by Lemma 2.3, there exists Ci and a £ R, such that 


E{\x{t)\^) = Y{t) < CiY{0)C = CiE{\xo\^)e, 


(3.12) 


where a satishes 2a + 6 + 2+ + (6 + 2d^)q°‘ = 0. According to a + 6 + + + < 0, we can know 

a < 0, the theorem is proved. 

□ 


Theorem 3.2. Suppose that x{t) is the solution of if the coefficients /, g satisfy assumption 

Id. 91 and 2a + 6 + 2c^ + (6 + 2df)/q < 0, then the zero solution of |d. 7| ) is almost surely polynomial 
stable. 

Proof. According to 2a + 6 + 2c^ + (6 + 2d^)/g < 0, we know that a + 6 + c^ + d^<0 and 
Li(|x(t)p) < CiE{\x{0)\‘^)C, where a < —1. 

By Ito formula and assumption 13.91 one can show that for any n — 1 < t < n. 


E{ sup \x{t)\'^) < E{\x{n — + {b + 2dJ)E / |x(gs)|^ds + 2c^£' / |x(s)pds 


n—l<t<n 


‘n—l 


' n—l 


+ 


2E{ sup f x'^{s) ■ g{s,x{s),x{qs))dB{s)). 

n—l<t<n Jn—l 


(3.13) 
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According to Burholder-Davis-Gundy inequations, it is easy to show that 


2E{ sup / x{sf-g{s,x{s),x{qs))dB{s)) 

n—l<t<n Jn—1 


rn 


< 8V2E{ / |a:(s)p|ff(s,a:(s),x(gs))pds)2 

J n—1 


f 1 

< 8\/2E( sup \x{t)\‘^ / |( 5 r(s,x(s),x(gs))pds )2 

n—l<t<n J n—1 

1 

<E{- sup |x(t)p + 64/ |sr(s,3:(s),x(gs))pds) 

^ n—l<t<n J n—1 

1 

<-E{ sup |x(t)p) + 64£'( / (c^lx(s)p + d^lx(gs)p)ds). 

2 n-l<t<n 


(3.14) 


' n—1 


Substituting (|3.14p into (|3.13|) . then 


E{ sup |x(t)p) <2C'i£'(|xoP)(n — 1)" + 132c^C'i£i(|xoP) / s“ds 


n—l<t<n 


'n—1 


pn 

+ {2b + l‘i2(f)CiE{\xo\^)q^ I s"ds. 

J n—1 

It is easy to know s“ds < (n — 1)“ max(l, 2“). So 

E{ sup |x(t)P) < (5(71 — 1)". 

n—l<t<n 

where C = [2 + 132c^ + (26 + 132d^)g“](i7ii?(|a;oP) ™ax(l, 2“). By Markov’s inequations, for any 
e > 0, it is not difficult to show 


P{ sup \x{t)f{n — iy 

n—l<t<n 


<- 


1 1 


1 


;E{ sup |x(t)P) 


<- 


^7 (n- 1)1+^ (71-1)“ Vi<t<n' 

‘ ' -c 


‘7 (ti — 1)1+'^ 

By using Borel-Cantelli lemma, the following limit can be achieved 

limsup sup |x(t)|^(7i — 1)“^“““^ = 0 a.s.. 

n^oo n—l<t<n 

Note that for any t > 0, there exists n{t) such that n{t) — 1 <t < n{t), and 

t 


lim 

t->oo n[t) — 1 


= 1 . 


Hence 


limsup \x{t)ft 

t—^OO 


2^—1—a—e 


< lim sup ( sup |x(s)|^(7i(t) — 1) ^ " ^)limsup 


t 


—l—a—£ 


t—^oo n{t)—l<s<n{t) 


t^oo in{t) - 1) 


—l—a—£' 


t 


—l—a—£ 


= lim sup ( sup |x(s)f(7i—1) “ ^)limsup 

t^oo n—l<s<n t—>-oo (^(v 1) 

= 0 . 


—j-) a.s. 

—l—a—e' 
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Due to the arbitrary of e, this can imply 


lim sup 
>-00 


iog\x{t)\ 

logt 


< 


1 + a 
2 


□ 


4 Consistence and convergence of the semi-implicit method 

In this section, we will employ the semi-implicit Euler methods to solve the equation (13.71) . 
We define a family of meshes with fixed step-size on the interval [0,T], i.e. 

T 

Tn = ■ ■ ,tN], fn = n/i, n = 0,1,2, • • • , V, /i = — < 1. (4.15) 

Since the points qt^ will probably not be included in T^v, so we need another non-uniform mesh 
which consists of all the points tn and qtn- Let 

Sn> = {0 = to = So, Si, S2, • • • , sn' = T}. (4.16) 

For any I, we have s/ = tn or si = qtm, where tn,tm £ Tjy. We can also display any si G S]\fi with 

tn< Si< tn+i by 

si=tn + Ch, Cs(0,1], t„,t„+iGTjv. 

In this paper, we denote by y{tn) the approximation of x{tn) at the point tn G Tjv, and y{qtn) 
the approximation of x{qtn) at the point qtn G Sn', then the semi-implicit Euler method is given 
by 


y{tn+l) = y{tn) + h[{l-e)f{tn, y{tn),y{qtn)) + Of{tn+l, y{tn+l), y{qtn+l))] +ff(in, y{tn) ,y{qtn))^Bn, 

(4.17) 

where y(to) = 2 : 0 , n = 0,1, 2, • • • , V — 1, = B{tn+i) — B{tn), 0 G [0,1]. Here we require y{tn) 

to be -measurable at the point n = 0,1, • • • ,N. 

We can also express (14.171) equivalently as 


y{tn+l) 


y{tn)+ / [{l-9)f{tn,y{tn),y{qtn)) + 0f{tn+l,y{tn+l),y{qtn+l))]<^t 

J tn 

rtn-\-l 

+ / g{tn,y{tn),y{qtn))dB{t). 

J tn 


(4.18) 


Note that we can’t express y{si) which equals to y{qtm) in (l4.18p . si G Sjy/, so we need a 
continuous extension that permits the evaluation of y{si)at any point si = tn + ^ ( G (0,1], 

so we define 


y(si) = y{tn) + [(1 - 0)f{tn, y{tn),y{qtn)) + 9fitn+l, y{tn+l),y{qtn+l))]dt 

+ / g{tn,y{tn),y{qtn))dB{t). 
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For any 0 given, tn & Tn, C ^ (0; 1]) the local truncation error of semi-implicit Euler method 
for (j.l.Tp can be denoted by 6h{tn,C)j 


Sh{tn, C) =x{tn + Ch) - {x{tn) + {1 - 0) 




f{tn,x(tn),x{qtn))dt 




+ ' 


f{tn+l,x{tn+l),x{qtn+l))dt 


(4.20) 


-h / g{tn,x{tn),x{qtn))dB{t)}. 

J tn 


Definition 4.1. (i) The semi-implicit Euler method is called to he consistent with order pi in 
average sense, if there exist constant C > 0 andpi, which are independent of step size h, such that 


max sup 0)1 < as h —)• 0. (4-21) 

o<’T’<^-iCe(o,i] 

(ii) The semi-implicit Euler method is called to he consistent with order p 2 in the sense of 
mean-square, if there exist constant C and p 2 , which are independent of step size h, such that 


max sup {E{\6h{tnX)\'^))^ <CW^ as h —)■ 0. (4.22) 

o<«<^^-iCe(o,i] 


For any 6 given, tn £ Tn, C £ (Ojl]) the global error of semi-implicit Euler method can be 
denoted by e(si) 


e{si) = e{tn + Qh) = x{tn + (h) - y{tn + Ch). 


(4.23) 


Definition 4.2. The semi-implicit Euler method is called to he convergent with order p, if there 
exist constant C and p, which are independent of step size h, such that 


max (E(|e(si)|^))2 < as/i —)• 0. (4.24) 

SiSSjv/ 


Lemma 4.1. Q/ Assume that there exists a positive constant K such that 

(i) (Lipschitz condition) Eor all t G [0,T], x\,X 2 ,yi,y 2 S 

\f{t,xi,yi) - f{t,X 2 ,y 2 )? V \g{t,xi,yi) - g{t,x 2 ,y 2 )? < - X 2 p + \yi - y 2 ?), 

(ii) (Linear growth condition) Eor all (t,x,y) G [0,T] x R'^ x , 

\f{t,x,y)(^ V \g{t,x,y)\‘^ < K{1 + \xf + |yp). 


Then there exists a unique solution x(t) to ((5.7|), and E{ sup \x{t)\‘^) < M. 

0<t<T 


Remark 4.1. Due to Lipschiz condition and /(t,0,0) = 0, g{t,0,0) = 0, it is not difficult to know 
\f{t,x,y)\^ < K{\x\'^ + lyp) and \g{t,x,y)\'^ < K{\x\'^ + \y\^). 


Theorem 4.2. Under the Lipschitz condition, the semi-implicit Euler method for equation (5.7) 
is consistent (i) with order 1.5 in average sense; (ii) with order 1 in mean-square sense. 
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Proof, (i) For the equation (13.7p and the semi-implicit method (14.171) . the local truncation 
error takes the special form: 


ShitnX) =(1-0) 




- f{tn,x{tn),x{qtn)))dt 






+ 


{fit,x{t),x{qt)) - f{tn+l,x(tn+l),x{qtn+l)))dt 


(4.25) 




{g{t, x{t),x{qt)) - g{tn, x{tn), x{qtn)))dB{t), 


for n = 0,1,2, • • • , A^. we will frequently make use of Hoder inequality in the next content. Note 
that F'dxl) < (Fidxp))^, so taking expectation and absolute both sides of the equation above, we 
can estimate 

|^(4(in,C))l <(1 - 0) / E{\f{t,x(t),x{qt)) - f{tn,x(tn),x{qtn))\)dt 

J t-n. 


(4.26) 


ptn-\-(^h 

+ 0 / E{\f{t,x{t),x{qt)) - f{tn+l,x{tn+l),x{qtn+l))\)dt 

•J tn 

, ftn + Ch 

<{l-e)K^ / [E{\x{t) - x{tn)\'^) + E{\x{qt) - x{qtn)\^)]^dt 

^ j-tn+Ch J 

-h OK^ / [E{\x{t) - x{tn+i)\^) + E{\x{qt) - x{qtn+i)\^)\'^dt. 

J tn 

By Lemma 14.171 h <1 and the integral we can obtain 

E(|x(t) - x{tn)\^) < {2K{t - tn) + 2K)E{ f (|x(t)p + \x{qt)\^)dt) < 8KMh. 

J tn 

In the same way, we can compute E{\x{qt) — x{qtn)\‘^) < 4:KMq{l + q)h, E{\x{t) — x(t„+i)p) < 
8KMh and E{\x{qt) — x{qtn+i)\^) < 4iFMg(l -|- q)h. Thus 


rtn+Ch ^ 

\E{6hitn,C))\< Cih2dt = CihH<Cih 

J tn 


(4.27) 


where Ci = K{8M + AMq{l + g)) 2 . 

This implies 

3 

max sup \E{5h{tn-,0)\ < Cih^. 

0<tn<N-l (g(0,l] 

(ii) According to the definition of Sh{tn,C)j fhe following inequality holds. 


rtn H”C, il 

E{\5h{tn, C)P) <3(1 - 0)^£^(| / {f{t, x{t),x{qt)) - f{tn, x{tn),x{qtn)))dt\‘^) 

Jtn 

+ 3B‘^E{\ / {f{t,x{t),x{qt)) - f{tn+l,x{tn+l),x{qtn+l)))dt\‘^) 

J t-n. 


+ 3E{\ 


f'tn-\-C^ 


{g{t,x{t),x{qt)) - g{tn,x{tn),x{qtn)))dB{t)\‘^) 






<3if((l - efch + 1)E{ / (|x(t) - x(t„)p + \x{qt) - x{qtn)\^)dt) 

J tn 

P “t“C^ 

+ 3K9‘^ChE{ / (lx(t) - x(tn+i)P + - x{qtn+i)\‘^)dt) 

J tji 

<C2h\ 

where C 2 = 2AK‘^M{9‘^ — 9 + l)(g^ + q + 2). Let C 3 = \/C^, then 

max sup (E(|4(tn,C)l^))^ < 

0<t„<Af-l (^g(o,l] 


(4.28) 


□ 


Theorem 4.3. Under Lipschitz condition, the semi-implicit Euler method for problem is 

convergent with order 0.5. 

Proof. For any si = tn + fh ^ Sn' , set 
Vh{tn, C) =(1 -G) x{tn),x{qtn)) “ f{tn, y{tn), y{qtn)))dt 

J tn 




+ ' 


{f{tn+l,x(tn+l),x{qtn+l)) “ f {tn+ 1 , y(tn+l), y{qtn+l)))dt 


(4.29) 


+ / {9{tn,x{tn),x{qtn)) - g{tn,y{tn),y{qtn)))dB{t). 

J tn 


Then 


e{si) = x(tn + Ch) - y{si) = e{tn) + dh{tn, C) + Vh{tn, 0- (4.30) 

Squaring both sides of the equation above, employing the conditional expectation with respect 
to the u-algebra Eq, and taking absolute values, we get 

^(|e(s/)| Vo) <S(|e(tn)| Vo) + E{\5h{tn, C)l Vo) + E{\vh{tn, C)l Vo) 

+ 2E{\e{tn)\^5h{tn.O\m + 2E{\e{tn)\-\MtnX)\\Eo) 

+ 2E{\5h{tnX)\ ■ Xh{tnX)\\Eo) a.s. 

=Ai + ^2 + ^3 + ^4 + A + ^6- 

Next we will estimate the six terms in ()4.28l) . For the term A 2 , by (I4.28P we have 
A 2 = E(lSh(tn, 01 Vo) = E(E(l6h(tn, 01 Vn)|-^o) < C2h^- 
For A^ in ()4.31l) . we obtain 


(4.31) 


2 I 3 <3(1 - 9fKChE{ 


f*^n+C^ 


{\x(tn) - y(tn)P + \x{qtn) “ y{qtn)\'^)dt\Eo) 


J tn 


+ 39^KChE{ / (lx(tn+l) - y(tn+l)P + \x{qtn+l) - y{qtn+i)\‘^dt\Eo) 

J tn 

+ 3KE{ / i\x(tn) - y(tn)P + \x{qtn) “ y{qtn)\‘^)dt\Eo) 

J tn 

=3KCh{{l - OfCh + l).F(|e(t„)| Vo) + 3KQh{{l - Offh + l).F(|e(gt„)| Vo) 
+ 30 V'h2.E(|e(t„+i)| Vo) + 30‘^KQ^h^E{\e{qtn+i)\‘^\EX. 
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We estimate A 4 , 

A 4 < 2(E (I E( 6 h(tn, 0)1^ 1 ^ 0 ))^ • (^(|e(t„)|Vo))^ < Cfh^ + hE{\e{tn)\^\Eo). 


In the same way, we can see 

^ <(2 - 9)KkhE{\e{tn)\^\Eo) + (1 - e)KkhE{\e{qtn)f\Eo) 

+ eKkhE{\e{tn+i)\^\To) + eKkhE{\e{qtn+i)\^\To), 

and 

Ae <2{E{\uh{tn,C)\^\m'^ • (^(|4(in,C)lVo))^ 

<C 2 h^ + 3KCh{{l - efCh + l)^(|e(t„)| Vo) 

+ W^KC^h^E{\e{tn+i)\^m 

+ mch{{l - efCh + l)E(|e(gtn)| Vo) 

+ 30V'h2^(|e((7tn+i)lVo). 

Combining these results, we can compute 

E{\e{tn + VI Vo) <(1 + QKCh{{l - efCh + l) + h + {2- 0)Kk/i)^(|e(tn)| Vo) 
+ {GO^KC^h^ + 9KHh)E{\eitn+i)\^\J^o) 

+ {6KCh{{l - 9)\h + !) + (!- 9)K"2Ch)E{\e{qtn)\^\Fo) 

+ {G9^KC^h^ + 0i^K/i)^(|eVn+i)l Vo) + ( 2 C 2 + Cf)h\ 

Set Rq = 0, Rn = max sup £^(|e(ti + Ch)p|-Fo)) then 

0 <i<n^g( 0 ,l] 

-E'(|e(tn)|^|-^o) < Rn, -E(|e(q'tn)Pl-^o) < Ru¬ 


in (|4.32l) . we need to calculate E{\e{qtn+i)\^\Fo) , which depends on either tn < qtn+i < tn+i or 
Q^n+l ^ tn- 

Case 1: If < Qtn+i < tn+i, then £l(|e((;t„+i)p|J^o) < Rn+i- According to ()4.32p . we can see 


^(|eV + Vl^l-^o) <(1 + 12(1 - 9fKh^ + l2Kh + h + (3 - 29)K^h)Rn 
+ {129^Kh‘^ + 29K^h)Rn+i- 


(4.33) 


So 


Rn+i = max sup E{\e{ti + Ch)\‘^\Fo) 

0 <»<»^+i ce(o,i] 

<(1 + 12(1 - 9fKh^ + l 2 Kh + h + (3 - 29)K^h)Rn 
+ {129^Kh^ + 29K^h)Rn+i + { 2 C 2 + Cf)h^. 


There is /iq = 2 ^ such that 1 — 129^Kh? — 29K‘2h > 0 when 0 < h < ho. Therefore 


^ ^ , 1 +12(1 - 0)V+ i2a: + 3aV 120 V, „ 

■^n+1 ^ 1 )^n “1“ 

1 - 129^Kh^ - 29K2h 


2 C 2 + Cl 


1 - 129‘^KK^ - 29K2h 


-h^ 


10 









Case 2: If qtn+i < tn, then E{\e{qtn+i)\‘^\To) < Rn- In the same way as case 1, we can get 


„ ^ i + i2{i-efK+ 12K+ m^+ i2e'^K.^ , 2 ^ 2 + <712 

Rn+l S (I + h-;;---i- )Rn + 


1 - 66»2Kh2 - eK2h 


- —h\ 

1 - 66»2Kh2 - OK^h 


when 0 < h < /ii = 2 . Now take 0 < L < 1, which is independent of h, such that 129“^Kh? + 

26K2h < L, and set 

1 + 12(1 - e)^K + 12i^ + 3Kh + ue^K ^ + 


l-L 


1-L 


Then combining case 1 and case 2, the Rn+i satisfies 

Rn+i <(1 + hM{9))Rn + C{e)h^ < (1 + hM{e))Rn-i + (1 + hM{e))C{e)h^ + C{9)h^ 

< • • • < (1 + hMi9))^+^Ro + c{e)h^ ^{1 + hM{e)y < 

1=0 ^ ' 

The expression above indicates that 

TM{0) _ 1 

^(|e(50l Vo) < ^n+i < -J^^C{9)h, 

for any si = tn + (h & Sn/ holds, where tn € T^r, ^ £ (0,1]. By the definition of convergence, we 
can show 


max (^(|e(si)|Vo))2 < 

siGS^f 


IC(9)(e™(f>^ - 1) 1 
M(0) 


The theorem is proved. 


□ 


References 

[ 1 ] Appleby, J.A.D., Berkolaiko, G., Rodkina, A., Non-exponential stability and decay rates in 
nonlinear stochastic difference equations with unbounded noise. Stochastics: An International 
Journal of Probability and Stochastics Processes, 81(2009), 99-127. 

[2] Appleby, J.A.D., Buckwar, E., Sufficient conditions for polynomial asymptotic behaviour of 
the stochastic pantograph equation. Stochastic Anal, 2003. 

[3] Appleby, J.A.D., Mackey, D., Almost sure polynomial asymptotic stability of stochastic dif¬ 
ference equations, Journal of Mathematical Sciences, 149(2008), 1629-1647. 

[4] Appleby, J.A.D., Mackey, D., Polynomial Asymptotic Stability of Damped Stochastic Dif¬ 
ferential Equations, Electronic Journal Qualitative Therory of Differntial Equations. 2(2004), 
1-33. 


11 










[5] Baker, C.T.H., Buckwar, E., Continuous 0-Methods for the Stochastic Pantograph Equation, 
Electronic Transactions on Numerical Analysis, 11(2000), 131-151. 

[6] Carr, J., Dyson, J., The functional differential equation y'{x) = ay{\x) + by{x), Proc. Roy. 
Soc. Edinburgh Sect. A. 74(1974), 165-174. 

[7] Ean, Z.C., Liu, M.Z., Cao, W.R., Existence and uniqueness of the solutions and convergence 
of semi-implicit Euler methods for stochastic pantograph equations, J. Math. Anal. Appl. 325 
(2007) 1142-1159. 

[8] Ean, Z.C., Song, M.H., Liu, M.Z., The ath moment stability for the stochastic pantograph 
equation,Journal of Computational and Applied Mathematics, 233(2009), 109-120. 

[9] Liu, K., Mao, X.R., Large time decay behavior of dynamical equations with random pertur¬ 
bation features. Stochastic Analysis and Applications, 19(2001), 295-327. 

[10] Liu, M.Z., Yang, Z.W., G.D.Hu, Asymptotical stability of the numerical methods with the 
constant stepsize for the pantograph equation, BIT, 45(2005), 743-759. 

[11] Mao, W., Convergence analysis of semi-implicit Euler methods for solving stochastic equations 
with variable delays and random Jump magnitudes. Journal of Computational and Applied 
Mathematics, 235(2011), 2569-2580. 

[12] Tsoi, A.H., Zhang, B., Weak exponential stability of stochastic differential equations. Stochas¬ 
tic Analysis and Applications, 15(1997), 643-649. 


12 


